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Summary 

A  mathematical  description  of  coupled,  surface-wave  propagation  in  elastic 
waveguides  whose  properties  change  slowly  in  the  lateral  and/or  transverse  di¬ 
rections  has  been  constructed.  This  construction  includes  curved  waveguides 
(shells)  whose  radius  of  curvature  is  large  with  respect  to  thickness.  The  work 
was  undertaken  to  explore  a  nondestructive  way  to  interrogate  an  inaccessible 
side  of  a  structure  by  launching  a  surface  wave  on  the  other,  accessible  side. 
For  flat  waveguides  (plates)  whose  thickness  varies  slowly  in  the  propagation 
direction  the  surface  waves  couple  very  efficiently  provided  the  thickness  is  on 
the  order  of  a  wavelength,  but  as  the  waveguide  grows  in  thickness  the  coupling 
is  quickly  lost.  The  coupling  in  a  curved  waveguide  occurs  provided  the  ratio  of 
thickness  to  radius  of  curvature  is  of  the  order  of  1/20  and  the  thickness  is  on 
the  order  of  a  wavelength.  The  modifications  caused  by  the  varying  geometry 
enter  the  form  of  the  guided  wave  as  slowly  varying  complex  amplitudes  and  do 
not  substantially  affect  the  dispersion  relation,  which  remains  very  close  to  that 
of  a  flat  plate.  Measurements  to  verify  the  existence  of  the  coupling  between 
surfaces  by  a  surface  wave  were  made  in  almost  flat,  metal  plates.  The  coupling 
was  found  by  Fourier  transforming  the  measurements  of  surface  displacement 
to  find  the  dominant  frequency,  which  represented  the  inverse  of  the  coupling 
distance.  Further  development  is  needed  to  transform  this  work  into  a  reliable 
nondestructive  evaluation  technique. 

Objectives  of  the  work 

1.  Develop  an  asymptotic  theory  for  time- harmonic  disturbances  that  de¬ 
scribes  the  propagation  of  the  two  lowest  waveguide  modes  in  a  curved 
elastic  structure  whose  radii  of  curvature  are  large  and  vary  slowly,  and 
explore  how  similar  they  are  to  those  of  a  flat  plate. 

2.  Explore  under  what  conditions  the  propagation  of  the  two  lowest  waveg¬ 
uide  modes  of  a  structure  can  be  described  as  the  propagation  of  a  Rayleigh 
surface  wave  along  one  or  the  other  of  the  two  bounding  surfaces  and  ex¬ 
plore  the  coupling  that  takes  place  between  the  surfaces. 
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3.  Show  that  using  the  coupled  elastic  surface  waves  is  one  way  to  inspect 
the  inner,  inaccessible  surface  of  a  structure  for  damage  from  the  outer 
accessible  surface. 

Relevance  to  the  Air  Force 

The  increasingly  longer  service  life  of  Air  Force  and  civilian  structures  means 
that  they  must  be  continuously  monitored  for  damage  to  ensure  their  reliability. 
A  surface  wave,  in  contrast  to  a  bulk  wave,  coupling  to  the  interior  surface  of 
a  plate  or  shell  would  strike  a  small  surface-breaking  crack  broadside  or  would 
remain  in  continuous  contact  with  a  patch  of  corrosion  and  hence  be  strongly 
perturbed.  This  would  indicate  the  presence  of  such  damage  at  an  early  stage 
before  it  became  severe,  and  hence  allow  the  structure  to  be  modified  or  removed 
from  service  before  it  failed. 

Outcomes 

Variable  thickness  waveguide .  At  high  frequencies  it  is  almost  impossible  to  ex¬ 
cite  either  of  the  two  lowest  Rayleigh-Lamb  modes  in  an  elastic  plate  separately. 
Instead  both  modes  are  excited  and,  because  the  lowest  is  antisymmetric  and 
the  next  lowest  symmetric,  the  two  propagate  together  and  mimic  a  Rayleigh 
surface  wave  that  couples  back  and  forth  between  one  surface  and  the  other.  We 
measured  this  phenomena  to  check  that  it  did  indeed  exist.  This  measuremnt 
work  is  reported  in  Ti  et  al  (1997).  We  were  successful  in  finding  the  phenomena 
though  its  measuremnt  was  not  as  direct  as  we  had  at  first  thought  it  would  be. 
A  measurement  of  the  particle  displacement  suggested  the  coupling  phenomena 
but  was  not  conclusive.  However,  Fourier  trasforming  the  record  of  the  particle 
displacement  clearly  indicated  the  coupling  from  one  surface  to  the  other.  That 
is,  there  appeared  a  dominant  frequency  in  the  measured  data  that  indicated  a 
steady  rise  and  fall  in  the  particle  displacement  and  thus  indicated  the  coupling. 
This  dominant  frequency  agreed  with  that  predicted  by  the  theory,  though  it 
did  not  appear  as  a  sudden  sharp  peak  but  as  a  smooth  hump  indicating  that 
the  coupling  was  not  as  clean  as  the  theory  would  suggest.  Nevertheless,  this 
suggests  how  a  nondestructive  testing  measurement  should  be  made:  one  should 
look  at  the  measurements  in  the  frequency  domain ,  rather  than  in  the  spatial 
one . 

We  undertook  to  calculate  this  coupling  phenomena  in  plates  with  variable 
thickness.  Figure  1  shows  the  problem  of  interest,  while  the  method  of  approach 
is  outlined  in  Figure  2.  The  asymptotic  expansions  used  have  the  form 

(1) 

Ti  =  e»W*'£,u> 
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Figure  1:  A  wedge-transducer  excites  the  two  lowest  Rayleigh-Lamb  modes  in 
an  elastic  plate  whose  thickness  varies  slowly  with  respect  to  wavelength. 
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Figure  2:  A  schematic  outline  of  the  mathematical  approach  taken.  The  term 
exp  means  exponential  and  contains  the  rapidly  oscillating  part  of  the  wave. 
The  modulated  amplitude  is  expanded  in  Rayleigh-Lamb  modes  at  each  point 
along  the  waveguide.  The  modes  are  those  corresponding  to  the  thickness  at 
that  point. 
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Figure  3:  The  normal  particle  displacement  is  shown  for  a  sinusoidally  varying 
waveguide.  The  magnitude  of  the  displacement  is  exaggerated.  The  frequency 
of  the  thickness  variation  is  very  much  smaller  than  the  frequency  of  the  elastic 
waves. 


Setting  the  coefficient  of  each  power  of  5  to  zero  gives  the  recursive  system 


This  formulation  of  the  eigenvalue  problem  for  a  Rayleigh-Lamb  modes  is  not 
widely  known  or  used.  The  t*  are  the  components  of  traction  and  the  Ui  those 
of  displacement.  The  a,  6  and  c  are  combinations  of  elastic  constants,  z  is 
the  transverse  coordinate;  x  the  lateral  one.  This  formulation  can  be  extended 
to  solve  a  very  broad  class  of  elastic  waveguide  problems.  The  details  of  this 
calculation  are  given  in  the  two  papers  Folguera  and  Harris  (1999a,b). 
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provided  the  caustic  remains  outside  the  upper  surface 
Rayleigh-like  surface  waves  couple  from  surface  to  surface 


Figure  4:  The  modes  of  a  cylinder  are  used  to  approximate  those  of  a  shell  with 
a  more  general  curvature. 

Figure  3  shows  our  results  for  a  plate  whose  surfaces  vary  sinusoidally.  The 
figure  exagerates  the  particle  displacement,  which  is  3£(u3),  the  real  part  of  the 
particle  displaceent  normal  to  the  surface. 

Curved  waveguides.  I  have  been  calculating  the  propagation  of  surface 
waves  on  an  arbitrarily  curved  shell-like  structure  by  trying  to  imitate  the  cal¬ 
culation  just  described.  My  first  attack  on  the  inplane  elastic  problem  for  the 
curved  structure  proved  unsatisfactory.  I  backed  away,  attacked  an  analogous 
scalar  problem,  solved  that  quite  satisfactorily  and  returned  to  the  inplane  elas¬ 
tic  problem.  This  new  calculation  is  now  almost  complete,  though  no  numerical 
calculations  of  results  have  been  done.  Figure  4  represents  geometrically  the 
general  scheme  of  the  approximation,  while  Fig.  5  tries  to  capture  the  mathe¬ 
matical  construction.  There  are  two  small  parameters  in  the  problem,  namely, 
*,  the  curvature  scaled  with  the  wavelength  and  S  the  rate  of  change  of  the 
curvature  over  a  wavelength.  Approximating  in  5  leads  to  a  recursive  sytem 
of  equations  similar  to  (2)  and  an  eigenvalue  problem  for  the  thickness  similar 
to  (3).  The  eigenvalue  problem  requires  finding  the  transverse  eigenmodes  for 
the  corresponding  cylindrical  structure  (Fig.  4),  which  leads  to  combinations  of 
Hankel  functions.  Rather  than  proceed  in  this  way,  the  eigenvalue  problem  is 
attacked  directly  using  an  expansion  in  k.  It  emerges  from  the  analysis  that  the 
product  of  thickness  times  the  curvature  htt  becomes  an  important  parameter. 
Two  limiting  cases  of  the  eigenvalue  problem  are  being  worked  out,  namely, 
hn  «  1  and  hn  »  1 

The  principal  conclusions  that  I  have  reached  so  far  are  the  following. 

1.  The  presence  of  any  curvature  affects  the  mode  shape  and  modulated  am¬ 
plitude  more  than  it  does  the  dispersion  relation.  The  dispersion  relation 
is  almost  identical  to  that  of  a  flat  plate  for  large  (that  is,  several  shear 
wavelengths)  radii  of  curvature. 
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Figure  5:  A  schematic  representation  of  the  structure  of  the  expressions.  The 
two  small  parameters  are  8,  the  rate  of  change  of  the  curvature  with  wavelength, 
and  k,  the  curvature  relative  to  the  reciprocal  wavelength.  Note  that  in  contrast 
to  Fig. 2  an  additional  expansion  is  neededto  solve  the  eigenvalue  problem. 

2.  In  any  curved  structure,  the  guided  waves  form  a  caustic  (each  guided 
mode  can  be  represented  by  a  set  or  pencil  of  rays;  the  rays  form  an 
envelope  called  a  caustic).  In  the  scalar  problem  the  caustic  can  never  be 
forced  to  lie  beyond  the  outer  concave  surface1,  though  it  can  be  forced 
to  lie  outside  the  inner  convex  surface  . 

3.  For  a  certain  range  of  the  parameters  the  caustic,  for  the  inplane  elastic 
problem,  can  be  forced  out  of  the  guide  into  a  region  beyond  the  outer 
concave  surface,  as  shown  in  Fig.  4.  In  this  case  coupling  of  Rayleigh 
surfaces  waves  between  the  surfaces  occurs  more  or  less  as  it  does  in  the 
flat  plate,  provided  Hk  «  0.1  . 

In  short ,  there  is  a  window,  largely  defined  by  the  frequency ,  shear-wave  speed, 
curvature  and  thickness  in  which  the  coupling  of  a  surface  wave  from  the  inner 
to  the  outer  curved  surface  occurs.  The  main  task  remaining  is  to  define  that 
window . 

Transitions  to  industry  /  Future 

I  have  become  involved  with  the  Technical  Center,  the  Advanced  Materials 
Group,  Caterpillar,  Inc.  The  mathematical  approach  developed  for  the  curved 
waveguide  work  is  now  being  used  to  study  propagation  in  thermal  barrier 
coatings.  These  coatings  are  used  in  engine  components  to  allow  the  engine  to 

1  Concave  and  convex  are  determined  by  imagining  yourself  within  the  guide 
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be  operated  at  higher  temperatures.  The  guided  waves  will  be  used  to  assess  the 
consistency  in  growing  these  coatings  and  to  monitor  the  their  integrity  when 
the  engine  is  in  service.  In  the  future  (next  two  years  or  so)  I  intend  to  continue 
to  study  propagation  in  coatings.  Please  do  not  hesitate  to  contact  me  if  you 
wish  to  learn  more  about  this  work. 

Administration  of  the  grant 

Starting  in  the  fall  of  1997,  Dr.  Alejandra  Folguera  was  supported  fully  on 
the  present  grant  as  a  post-doctoral  associate.  Starting  in  September  1998, 
Dr.  Folguera  was  supported  one-third  time  on  this  grant,  while  teaching  in  the 
Mathematics  Department,  UIUC,  for  the  other  two- thirds.  In  May  1999  she 
choose  to  end  her  appointment.  With  what  funds  remained  I  supported  myself 
over  the  summer  of  1999.  Two  publications  acknowledging  this  grant  have 
appeared,  namely  Folguera  and  Harris  ( 1999a, b).  A  further  paper  describing 
the  curved  surface  work  will  be  submitted  for  publication  and  a  shorter  paper 
descibing  this  work  will  appear  in  a  conference  proceeding. 
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